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According to recent Quantum Monte Carlo simulations the 
small polaron theory is practically exact in a wide range of 
the long-range (Frohlich) electron-phonon coupling and adi- 
abatic ratio. We apply the Lang-Firsov transformation to 
convert the strong-coupling term in the Hamiltonian into the 
form of an effective hopping integral and derive the single- 
particle Green's function describing propagation of the small 
Frohlich polaron. One and two dimensional spectral func- 
tions are studied by expanding the Green's function pertur- 
batively. Numerical calculations of the spectral functions are 
produced. Remarkably, the coherent spectral weight (Z) and 
effective mass {Z') renormalisation exponents are found to be 
different with Z' » Z, which can explain a small coherent 
spectral weight and a relatively moderate mass enhancement 
in oxides. 

PACS numbers:74. 20. Mn,74.20.-z, 74.25. Jb 

1. Introduction 

The problem of a fermion on a lattice coupled with the 
bosonic field of lattice vibrations has an exact solution in 
terms of the coherent (Glauber) states in the extreme 
strong-coupling limit, X ^ oo for any type of electron- 
phonon interaction jicmserving the on-site occupation 
numbers of fermionsBH. For the intermediate coupling 
the 1/A perturbation diagrammatic technique has been 
developed both for a singleQ and multi-polaix)n|Siistemsa. 
The expansion parameter actually is 1/2zAtJ'EII3ij, so the 
analytical perturbation theory might have a wider region 
of applicability than one can expect from a naive vari- 
ational estimate {z is the lattice coordination number). 
However, it is not clear how fast the expansion converges. 
The exact numerical diagonalisation of vibrating clus- 
ters, variational calculationsLTtZL dynamical mean- field 
approach in infinile jdimensionstS, and Quantum-Monte- 
Carlo simulationsEHrtj revealed that the ground state en- 
ergy (the polaron binding energy Ep) is not very sensitive 
to the parameters. On the contrary, the effective mass, 
the bandwidth and the shape of polaron density of states 
strongly depend on polaron size and adiabatic ratio in 
case of a short-range (Holstein) interaction. In particu- 
lar, numerical diagonalisation of the two-site-one-electron 
Holstein model in the adiabatic ojq / 1 < 1 a,s well as in the 
nonadiabatic uo/t > 1 regimes shows that perturbation 
theory is almost exact in the nonadiabatic regime for 
all values of the coupling constant. However, there is no 
agreement in the adiabatic region, where the first order 
perturbation expression overestimates the polaron mass 
by a few orders of magnitude in the intermediate coupling 
regimetl. Here ojq is the characteristic phonon frequency 
and t is the nearest-neighbour hopping integral so that 



the dimensionless coupling constant is A = Ep/{zt). A 
much lower effective mass of the adiabatic small polaron 
in the intermediate coupling regime compared with that 
estimated by first order perturbation theory is a result 
of poor convergence of the perturbation expansion owiug 
to the appearance of the familiar double- well potentials 
in the adiabatic limit. The tunnelling probability is ex- 
tremely sensitive to the shapCi-ef this potential. 

It has also been understoodEilptliat the range of appli- 
cability of the analytical theorydu strongly depends on 
the radius of interaction. While the analytical approach 
is applicable only ii ujq > t for short-range interaction, 
the theory appears almost exact in a substantially wider 
region of the parameters for a long-range (Frohlich) in- 
teraction. The exact effective mass of both small and 
large Frohlich polarons, calculated with the continuous- 
time path-integral Quantum Monte Carlo (QMC)_algo- 
rithm, to* (A) is well fitted by a single exponentcj'tJ. As 
an example, e"'^'^'^ for ujq = t and e^ '^^ for loq = 0.5 1 
describe to* (A) in a one-dimensional lattice. The numer- 
ical exponents are remarkably close to those obtained 
from the Lang-Firsov transformation, e°-^*'^ and e^'^^'*', 
respectively. Hence, in the case of the Frohlich interac- 
tion the transformation is perfectly accurate already in 
the first order of 1/A expansion even in the adiabatic 
regime, ujo/t < 1 for any coupling strength. 

In this paper we use this result to calculate the single- 
particle Green's function of a fermion on a lattice coupled 
with the bosonic field via the long-range Frohlich inter- 
action. 

2. Green's functions of the small Prohlich polaron 

The classical approach to the small polaron problem is 
based on the canonical displacement (Lang-Firspv) trans- 
formation of the electron-phonon HamiltonianB allowing 
for the summation of all diagrams including the vertex 
corrections, 

H = ^t(m - n)Ss^s'4^j + ^ Wqfii [ui(q)dq + H.c.] 

+ ^C.q(dtdq + l/2) (1) 

q 

with the bare hopping integral i(m) and the matrix ele- 
ment of the electron-phonon interaction 

^..(q) = ^^(q)e»q-. (2) 

Here i = (m, s), j = (n, s'), include site m and spin 
s, fii — clci, and Ci,dq are electron (hole) and phonon 
operators, respectively. 
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The canonical transformation e"^ diagonalising the 
Hamiltonian in the X = oo hmit is 



where 



S* = ^ hi [ui{q}dq - H.i 



(3) 
(4) 
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The electron operator transforms as 



Ci = CiCXp 



and the phonon one as: 

dq = dq + ^hiU*{q) 



(5) 



(6) 



It follows from Eq.(6) that the Lang-Firsov transforma- 
tion shifts ions to new equilibrium positions. In a more 
general sense it changes the boson vacuum. As a result, 

i.j i 

+ ^Wq(d]jdq + l/2) + -^Uyn,nj, (7) 

q i=i] 

where 

<T„ = i(m - n)Ss,s' exp ^^[uj(q) - Mj(q)]dq - H.c^ 

(8) 

is a renormalised hopping integral depending on the 
phonon variables, 



2N 



(9) 



the polaron level shift (binding energy), and 

l^|7(q)P^qCOs[q.(m-n)] (10) 



the attractive interaction of polarons owing to lattice de- 
formation. 

We consider a single-particle problem wheri-the last 
term in Eq.(7) is absent. The QMC resultEil allows 
us to average the first hopping term with respect to 
phonons no matter what the value of the Frohlich in- 
teraction is. Then we can apply the canonical trasforma- 
tion together with the averaged Hamiltonian to derive the 
single-particle Green's function (GF). The interaction is 
completely integrated out from the averaged Hamilto- 



H = Hp + Hph (11) 
where the 'free' polaron part is given by 

i/p = ^C(k)4ck, (12) 

k 

(from now we omit the spin index), and the free phonon 
part is 



Fpft =^c^q(4dq + l/2). 



(13) 



Here ^(k) — Z'i?(k) — /x is the renormalised polaron- 
band dispersion with the chemical poten|tial fi, which 
includes polaron binding energy {—EpjBi^ E{]i) = 
^jjj i(m)ea;p(— ik • m) is the bare disperison in a rigid 
lattice, and the mass-renormalisation exponent is 



with 



7/ _ Em i(m)e-f (™) exp(-ik • m) 
Em exp(-ik • m) 



g2(m)=^|7(q)ni-co5(q.a)]. 



(14) 



(15) 



Quite generally one finds Z' = exp{—^Ep/Lu), where the 
numerical coefficient 
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= J2 l7(q)Hl - cos(q • a)]/ ^ ^(q)^ (16) 



might be as small as 0.4E3 and even smaller in the 
cuprates with nearest neighbour oxygpni-oxygen distance 
less than the lattice constant, 7 ~ 0.2l§c3. 

Applying the Lang-Firsov canonical transformation 
the Fourier component of the retarded GF, 



^ "'0 



-if k-m— ojtl 



{0\co{t)cl{Q)\0)dt 



(17) 



is expressed as a convolution of the Fourier components 
of the coherent retarded polaron GF, Gp{n,t), and the 
multiphonon correlation function cr(n, t): 

^'^ „ J-00 



(18) 



Here 



Gp{nY,uj)^~i dte*"*(6|e'-f^''*coe-*^''*ct„|0), (19) 
Jo 



and 



CT(m,w)=/ dte''^\Q\e'"^'^'Xoe''"^'^'xm (20) 

JO 
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with — cxp ^X]q ''^i(q)'^q ^ -ff-c.^ , and |6) the ground 

state of H. IIP 
Straightforward calculationsoB yield 



1 



-m 



7V^u;-^(k')+i(5' 



(21) 



and 



CT(m, w) — iZ 



1 



(=0 



Here 



1^(^)1^ ^^P*^*^' 



Z = exp(-^|7(q)P). 



(22) 



(23) 



In the foUowing we consider dispersionless phonons, ojq = 
Wo, and the Frohhch interaction with 7(q) ~ 1/q. The 
convolution of Eq.(19) and Eq.(20) yields 



Gfl,(k,w) = ^G^)(k,w), 



(24) 



i=0 



where 

G^Hk,w) = 

V l7(qi)P X |7(q2)P X ... X |7(qOP 
^^^^^ (2iV)'n(w-wo?-C(k + qi + ...qO+*<5) ^ ^ 

with (5 = +0. 

Obviously, Eq.(22) is in the form of a perturbative mul- 
tiphonon expansion. A term with an index I corresponds 
to a transition from the initial state k of the polaron 
band to the final st^e k + qi + ...qi with the emission 
of I optical phononsEZl. 



3. Spectral functions of the small Frohlich polaron 

The Green's function of a polaronic carrier, Eq.(22) 
comprises two different contributions. The first (/ = 0) 
coherent k-dependent term arises from the polaron band 
tunneling. 



G 



(0) 



-m 



(26) 



The spectral weight of the coherent part is strongly (ex- 
ponentially) suppressed as Z — exp{—Ep/LUo) while the 
effective mass might only be slightly enhanced, — 
Z'E]^ — ^, because Z « Z' < \ in the case of the 
Frohlich interaction. 

The second incoherent phonon-assisted contribution 
with I > 1 describes the excitations accompanied by 
emission of phonons. We believe that this term is respon- 
sible for the background in optical conductivity and in 



photoemission spectra of cuprates and manganites. We 
notice that its spectral density spreads over a wide en- 
ergy range of about twice the polaron level shift Ep. On 
the contrary the coherent term shows an angular depen- 
dence in the energy window of the order of the polaron 
bandwidth, 2Z'zt. Interestingly, there is some k depen- 
dence of the incoherent background as well, if the ma- 
trix element ofj-eJectron-phonon interaction depends on 
q (see also Ref.Ea). To illustrate this point we calculate 
the single-phonon contribution {I = 1) to the spectral 
function 



A(k,w) = -(l/7r)3G;^(k,w) 



^«(k, 



(27) 



The coherent part {I — 0) of the spectral function is 
a (5-fupjction in agreement with the well-established fact 
( see0D and references therein) that small polarons ex- 
ist in the Bloch states at zero temperature no matter 
which values the parameters of the system take. The 
single-phonon contribution to the incoherent background 
is given by 



AW(k,w)~ / dqq-H{LJ-iJo~a^ + q), 
JBZ 



(28) 



where the interval of integration (the Brillouin zone, BZ) 
is determined by lattice constants a, b, c. We calculate 
this integral for one-dimensional (ID), ^(k) = 2tcos(fc:Ea) 
and two-dimensional (2D), ^(k) = 2icos{kxa) + 
2t' cos{kyb) polaron bands in the tight-binding approx- 
imation with the (renormalised) nearest neighbour hop- 
ping integrals, i = Z't, i' = Z't' . 

The ID single-phonon spectral function is reduced to 

A«(k,w)~ (1-^2^-1/2 
f-7r/3 

dz[/(z,w,A:,)-f /(z,w,-fc,)], (29) 

^0 

where 

/(z, w, fc,) = + [k, - cos-^Qf]-y^ 

tan"i( TT [2^ + (k.^ - cos-^Cjf]-^''^). (30) 

Here and further we take a = h = 1 and c = 3, and 
t' = t/4. Its spectral and angular dependence are shown 
in Fig.l. Apart from two nondispersive ID van-Hove 
singularities (vHs) at a) = (uj — uj^)/2t — ±1 there is 
an interesting dispersive peak aX uj ~ cosk^, which is 
due entirely to the long-range character of the Frohlich 
interaction. Indeed, approximating the Brillouin zone 
by a cylinder along x, one readely obtains a logarithmic 
singularity in the spectral function, 

AW (k, w) ^ (1 - w^)- V2 1, ( {k. -C0S-^^? + ql 

(31) 
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Here qo is the radius of the cyhnder. For uj = ujq + 
2tcoskx, the x-component of the phonon momentum is 
zero, and the singular matrix element squared 1/?^) 
integrated over qy and q^ yields a singularity. 

The 'long-range' dispersive features appear in the 2D 
single-phonon spectral fuction as well. This function is 
reduced to 

A(l)(k,C<j) - /''^"^dz[l- (W- icOSz)2]-l/2 

X [ (z, w, kx,ky) + (t){z, UJ, kx, -ky) 

+ (/)(z, w, -kx, ky) + (j){z, u, -kx, -ky)], (32) 

with 

(j){z, UJ, kx, ky) = [(.T - kx)"^ + fcy)^]~^^^ 

tan-i (^[(x - kx)^ + (z - ky)X^/^) (33) 

and X — cos^^(a) — ^cosz). The integration limits are 
defined as a{uj) = cos-i[4(c:; + 1)] if -5/4 < Q < -3/4, 
a(w) = if -3/4 < w < 5/4, and = cos-i[4(w - 1)] 
if 3/4 < ci) < 5/4, 6(cj) = TT if -5/4 < ci) < 3/4. The 
numerical results are shown in Fig. 2. 

We do not expect any dispersive peaks in the multi- 
phonon {I > 2) contributions to the spectral function 
because of the additional integrations compared with the 
single-phonon term. However, the nondispersive vHs re- 
mains for all ^ > 1. They can be washed out by the 
phonon frequency dispersion in real crystals. 

4. Conclusions 

We have calculated the Green's and spectral functions 
of a fermion on a lattice coupled with lattice vibrations 
via the long-range Frohlich interaction. In a sharp con- 
trast with the Holstein polaron the mass renormalisation 
exponent of Frohlich polaron {Z') differs from the renor- 
malisation of the coherent spectral weight (Z). On the 
one hand this important result tells us that small po- 
larons as well as intersite bipolarons are perfectly mobile 
and can account for the high-Tc values in cupratesE3. 
On the other hand the coherent spectral weight remains 
strongly suppressed in polaronic conductors, Eq.(24), be- 
cause Z might be less than Z' by one or even a few orders 
of magnitude. These unusual spectral features provide an 
explanation for the apparent discrepancy between a very 
small Drude weight and a relatively moderate mass en- 
hancement, TO* ~ 3mep=j_li)TOe (depending on doping) of 
carriers in manganitesEBEEl and cuprates. They also ex- 
plain why the.'p.ytended' vHs observed in angle-resolved 
photoemissior£lfE3 can be hardly seen in angle- averaged 
photocmission. Indeed, the integrated spectral weight of 
the incoherent background is proportional to {1 — Z). It 
turns out to be much larger than the coherent contribu- 
tion, proportional to Z << 1. Finally, the k dependent 
incoherent background. Fig. 1,2, might obscure the ex- 
perimental determination of a Fermi-level crossing. We 
believe that our results can provide a quantitative ap- 
proach to the experimental tunneling and photocmission 



spectra. In particular the coherent part of ARPES and 
tunneling in cuprates has been recently explained^. 

The authors greatly appreciate enlightening discus- 
sions with Chris Dent. C.S. has been supported in this 
work by a grant from the Royal Thai Government. 
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Figure Captures 

Fig.l. ID single-phonon contribution to the polaron 
spectral function. 

Fig. 2. 2D single-phonon contribution to the polaron 
spectral function along the T — Y direction. 
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